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Abstract 

The effect of polarization rotation on the performance of metal oxide semiconductor field-effect 
transistors was investigated with a Landau-Ginzburg-Devonshire theory based model. In this 
analytical model, depolarization field, polarization rotations and the electrostatic properties of the 
doped silicon substrate are considered to illustrate the size effect of ferroelectric oxides and the 
stability of polarization in each direction. Based on this model, we demonstrate that MOSFET 
operation could be achieved by polarization reorientation with a low operating voltage, if the 
thickness of ferroelectric oxide is properly selected. Polarization reorientation can boost the surface 
potential of the silicon substrate, leading to a subthreshold swing S lower than 60 mV/decade. 
We believe that this model could provide guidance in designing electronic logic devices with low 
operating voltages and low active energy consumption. 
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I. INTRODUCTION 


Ferroelectric oxides are a promising class of materials for application in electronic devices 


due to their intrinsic spontaneous electric polarization, which can not only contro 


the con- 


Uii. 


By 


ductance of the channel, but also can be reoriented by an external electric held 
modulating the polarization of ferroelectric oxides, programmable binary logic devices can be 
achieved, and the fast reorientation of polarization enables fast switching and lower-power 


, 0 - 


operation of the metal-oxide-semiconductor held-effect transistor (MOSFET) 

Here, we aim to provide guidance about designing a programmable fast switching MOS¬ 
FET with proper size, by considering factors which were rarely included in previous mod¬ 
eling, but may strongly affect the polarization reorientation and size effect of ferroelectric 
oxides. Many analytical models based on the Landau-Ginzburg-Devonshire (LGD) the¬ 
ory have been proposed previously [?, 0 to simulate the electrical behaviors of MOSFETs. 
These models provide an insightful understanding about the mechanisms of ferroelectric 
oxide based MOSFETs and guide the fabrication of novel devices. However, there are still 
several vital factors beyond the scope of previous models. First, the effect of the polarization 
distribution in three dimensions (3D) on the channel current-gate voltage relationship of a 
MOSFET was rarely considered, even though there were a lot of studies about polarization 


in 3D and its response to electric helds in different orientations jd, 9-|ll|. For simplicity, the 


polarization of the ferroelectric oxide in a MOSFET is usually treated in one dimension. It 
is true that the channel conductance is mainly modulated by the out-of-plane polarization, 
but it is also important to note that the polarization components in all three dimensions 
are coupled together, and the in-plane polarization strongly influences the electric suscep¬ 
tibility out of plane. The second factor is the electrostatic properties of the channel and 
gate electrode. It is widely known that the distribution of charge in electrodes, which is 
parameterized as screening length, determines the strength of depolarization held, which 
affects the magnitude of the spontaneous ferroelectric polarization 12l-ll5|. 

In this letter, we propose a LGD theory based single crystal model with detailed analysis 
of these factors, in order to provide strategies for designing a low operation voltage ferroelec¬ 


tric held-effect transistor [16 


, [l7| . A previous study argues that a subthreshold swing lower 


than 60 mV/decade can be achieved by the negative capacitance ehect 


17| , but there is also 


debate that direct current negative capacitance is not possible, due to Gibbs free energy con- 
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siderations 


18| . In our model, the polarization dynamics obeys the Landau-Khalatnikov 


equation and is always minimizing the Gibbs free energy under a unidirectional gate voltage. 
We demonstrate that fast switching (subthreshold swing lower than 60 mV/decade) can be 
achieved by a proper design of the ferroelectric oxide size. The mechanism is that during 
the process of polarization reorientation, the tendency to possess spontaneous polarization 
in a new direction boosts the screening charge accumulation and channel current increases, 
leading to a low subthreshold swing. 


II. MODEL APPROACH 


The LGD model is a phenomenological theory which describes the electrical properties 
of ferroelectric oxides. In this model, the thermodynamic potential (Gibbs free energy 
Go) of a single crystal ferroelectric oxide is given as a function of polarization in three 


directions 


19 


20 |. 


Go = ai (Pj + P; + P;) + an {P^ + P„" + P^) 

+ an (P|P; + PlPl + P^Pj) + am (P* + P® + P®) ^ 

+ a.io K (Pi + Pi) + Pi (Pi + Pi) + Pi (Pi + Pi)] 

+ aniPlPlPl 

Taking external electric held and internal depolarization held into consideration, electrostatic 
terms should be added as 


G — Gq — E^Px — EyPy — EzPz 


( 2 ) 


FIG. 1. shows the schematic of the MOSFET that we study. The insulator between 
the gate electrode and doped silicon substrate is a ferroelectric oxide jl^ . 1^ . The 2 ; axis 
is normal to the ferroelectric oxide/silicon interface. In the x and n directions, there is no 


external voltage, and short circuit conditions are applied 


[3,y, 


23|. For the case that a 


gate voltage Vg is imposed on the MOSFET, we have the following equations: 


214x + Vox = 0 

“^Vey + Voy = 0 ^ 

Vez + Voz + (fs = Vg ^ 


( 3 ) 
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p-type Silicon 



FIG. 1. The schematic of a MOSFET with ferroelectric oxide as the insulator between the gate 
electrode and silicon substrate. The pink rectangles represent insulator layers isolating the side 
electrodes and source, drain terminals, (a) No gate voltage is applied iVg = 0) and the polarization 
is in plane. No carriers or current is in the channel; (b) Gate voltage is applied {Vg > 0) and 
polarization is out of plane. Carriers are induced by the polarization, and drain-source current 
flows. 


Vex,ey,ez and Vox,oy,oz are the voltage drop across the electrode and the ferroelectric oxide in 
the X, y and z directions, ips is the snrface potential of the silicon snbstrate, and it can also 
be viewed as the voltage drop in the doped silicon snbstrate. The flat band potential V/b, 
which resnlts from the alignment of the Fermi levels of the gate electrode, oxide, and silicon 
snbstrate is inclnded in Vg. The electric held E is determined by both the external applied 
voltage and the electrostatic properties of the ferroelectric oxide and electrodes 24l-l26l|. It 
is widely accepted that the charge density in noble metal electrodes follows the Thomas- 
Fermi distribntion, and this distribntion canses a voltage drop across the electrodes. The 
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following derivation calculating this potential drop follows the main idea in Ref. 22|], but 
is re-interpreted. Taking the 2 ; direction as an example and the electrostatic properties in 
the X, y directions following similar rules, the relationship between the electric held and the 
charge density takes the form 


dE{z) n {z) - no 

= -Q{z) = -q- 


dz 




(4) 


E (z), Q (z) and n [z) are the electric held, the charge density and the electron density in 
electrodes at the position z. no is the average electron density in a neutral electrode, q is 
the electronic charge. £0 and Eg are the electric permittivities of the vacuum and electrode. 
Meanwhile, the potential drop V (z) is expressed as 


dV {z) 
dz 


dV (z) dn (z) 
dn {z) dz 


-E{z) 


(5) 


The electrons in metal electrodes are treated as a free Fermi gas, so the local potential and 
the electron density are related as 


V = (STT^n) 


dViz) 


2m 


(Svr^) ^ n{z) 3 


( 6 ) 


(7) 


dn (z) 3m 

h is the reduced Planck constant and m is the electronic mass. By combining equations (5) 
and (6), we have 

dn {z) 


dz 


2 I 

-(3vr^) ® n {z)~^ 

3m ^ ^ ^ 


-1 


Eiz) 


Taking the derivative of equation (4), we have 

d^E (z) q dn (z) 3mq 


dz‘^ 


EoEg dz EoEgK^ 


(dvr^) ^n{z)^E{z) 


( 8 ) 


(9) 


The characteristic length A^, which is also called screening length and determines the dis¬ 
persion of electrons in electrodes, is dehned as 


A 


2 

2 


3mg 

EoEgh? 





Therefore, equation (9) is rewritten as 


d^E (z) 
dz‘^ 


1 


E{z) 


( 10 ) 


( 11 ) 
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The boundary conditions are 


^ E (0) = 

£0^e 

( 12 ) 

E (— 00 ) = 0 

Qz is the screening charge density at the ferroelectric/oxide electrode interface, which is also 
Q (z = 0). Thus, the electric held and potential drop through one electrode are 


B (z) = 


rO 


Ve= / E{z)dz = 


9Ee-/^^dz = ^ 


(13) 


(14) 


— 00 J — 00 eQEf. e^Ef. 

The heterostructure of electrode, ferroelectric oxide, and silicon substrate can be regarded as 
a capacitor, with equal charge densities at each interface. However, the charge distribution 
in the doped silicon substrate is quite different from that in metal. This is because electrons 
in the metal are treated as a free electron gas. This is the basic approximation of the 
Thomas-Fermi model. But doped silicon is a semiconductor, and the free carrier density is 

271 . |28| . The interface charge density-potential relationship in the 


local potential dependent 
silicon substrate is given by 

Qz = \/2EsikTNa 

d(p{z) , l2kTNn 


_2i£s q(ps 
e fcT + 


dz 


= E{z) = 


£Si 


_ QVW 
e kT -\- 


kT 

kT 




-1 + 




iV2 


^ _ gif {z) 
kT 


- 1 


(15) 


(16) 


(Ps is the surface potential of the silicon substrate, k is the Boltzmann constant. Other 
parameters are listed and described in TABLE I. 

From the analysis above, we see that the charge density decreases gradually away from 
the oxide in both the metal electrode and the doped silicon substrate, even though the 
analytical expressions and physical mechanisms which govern the charge distribution are 
different. As a result, there are voltage drops through each layer. These voltage drops could 
counteract or completely neutralize the applied gate voltage, exerting signihcant influence 
on the magnitudes of ferroelectric polarization and charge in the channel. Equation (15) 
demonstrates that there is a one-to-one correlation between the interface charge density Qz 
and the surface potential ps- Ts is a function of Qz'- 


Ts = f{Qz) 


(17) 
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The voltage drop across the ferroelectric oxide takes the form 

K. = E,-d,= 

^0 


( 18 ) 


dz is the thickness of the ferroelectric him and Pz is the polarization in the z direction. With 
the analysis above, equation set (3) is rewritten as 


Qx^x ^ Qa 

£o£e 

Qy\ 


Qz^z 


+ 


+ 


^0 

Qy ~ Ey 


Px 7 , Qx^a 

-dx + 


Qy^y 


^0 

Qz — Pz 


= 0 


= 0 


dy + 

dz + f {Qz) = Vg 


(19) 


For short-circuit conditions, in order to balance the potential drop in the electrodes, the 


sign of Vox,oy should be opposite to that of Vex,ey This indicates that surface charge density 
should be smaller than polarization, which means an incomplete screening of the polariza¬ 
tion charge. As a result, an electric held (depolarization held) is induced opposite to the 
polarization. The potential drop in the metal electrodes, which is proportional to screening 
length, is the origin of the incomplete polarization charge screening and the depolarization 
held which suppresses ferroelectricity. 


The energy surface versus polarization direction and magnitude can be plotted under 
the electrostatic restrictions expressed in equation (19). After acquiring the energy sur¬ 
face, polarization dynamics on the energy surface is simulated by the Landau-Khalatnikov 
equation j29l-l8l|. 


7- 


d'^ 


dt 


V^G' = 0 


( 20 ) 


7 is the polarization dynamic parameter. G is the thermodynamic potential dehned in 
equation (2) with the restriction shown in equation (19). The most stable polarization is 
the one which minimizes Gibbs free energy G. However, if the polarization is not in a 
local minimum, it cannot move to one instantaneously. The rate of return to a minimum 
is determined by many factors. For example, the resistance of the circuits ahects this rate, 
because polarization evolution must be accompanied by screening charge transmission. The 
polarization dynamic parameter 7 is related to the mobility of polarization, as can be 

regarded as the driving force of polarization and 7 ^ is the speed of polarization evolution. 
The applied time-varying gate voltage takes the form. 


Vg = Vo sin (cut) 


0 < t < 


TT 


U 


( 21 ) 
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Here, we do not mean that the applied gate voltage is oscillatory. Instead, we are simulating 
one on/off programmable cycle (O < t < ^), and the increase/decrease of the gate voltage 
takes the sine form. Equation fl20l) is rewritten as 


dP 

d {ojt) 


+ — 0 


( 22 ) 


7 o = cjy is the effective polarization dynamic parameter, and Qz can be calculated from 
Pz and the drain-source current Ids is obtained by the Pao-Sah double integral j^. 

'■V^s 


T - ^ 

^ds 




(23) 


where the parameters in this simulation are listed in TABLE I, and the function ^ {ip, V) is 
given in Ref. 


III. RESULTS AND ANALYSIS 

The ferroelectric oxide we choose is BaTiOs, which possesses a relatively large sponta¬ 
neous polarization (P* ^ 0.26 C/m^) at room temperature jl^ . 

In order to simulate the energy surface, we vary the surface potential ps and polarization 
Px- For each ps, charge density Qz and polarization Pz are determined uniquely by equations 
(15) and (19). At room temperature, the BaTiOs crystal has a tetragonal phase. The 
polarization orients either out of plane or in plane. We set the in-plane polarization direction 
as the X direction and Py = 0. Here, we should also note that we assume that the in¬ 
plane polarization has no effect on the channel. Therefore, it is not necessary that the 
source channel-drain-current flows along the x direction. Electric held E is obtained by 
the electrostatic restrictions in equation (19). Then energy surfaces describing Gibbs free 
energy G with respect to Px and Pz are calculated by formula (1) and (2). In FIG. 2, we have 
plotted two energy surfaces of BaTiOs with different thicknesses in the x and z directions 
on a p-type silicon substrate. 

From the graphs, we can see that for out-of-plane polarization, a negative orientation 
(pointing to gate electrode, with negative ends of oxide dipoles toward the channel) is more 
favorable when there is no applied voltage. This is because for a p-type silicon substrate, 
positive screening charge is more likely to accumulate at the interface, leading to the po¬ 
larization pointing away from the substrate/ferroelectric oxide interface. FIG. 3 shows the 





TABLE I. PARAMTERS INVOLVED IN THIS STUDY 



Description 

Value 

T 

Temperature 

298 K 

ai 

Coefficient in LCD theory^ 

-2.77x10^ m/F 

ail 

Coefficient in LCD theory^ 

-5.35x108 mVC^F 

012 

Coefficient in LCD theory^ 

3.23x10® mVC^F 

am 

Coefficient in LCD theory^ 

8 .00x10^ m^/C^F 

0112 

Coefficient in LCD theory^ 

4.47x10^ mVC^F 

0123 

Coefficient in LCD theory^ 

4.91 xlO^ mVC^F 

^x,y,z 

Screening lengths in noble metal^ 

0.04 nm 

£ 

Dielectric constant of noble metah 

2.0 

Na 

Substrate doping concentration*^ 

4 X 10^^ cm“8 

m 

Intrinsic carrier concentration*^ 

1.5 X 10^1 cm"® 

£Si 

Dielectric constant of silicon 

11.7 F/m 

/^eff 

Effective electron mobility 

3.0 X 10“^ m^/Vs 

w 

Width of the silicon channel 

4.0 X 10“^ m 

L 

Length of the silicon channel 

4.0 X 10“^ m 

dy 

Equals to dx 



Reference 

[R 

Reference 

26 

Reference 

28 


relationship of the surface potential and the interface charge density in the p-type silicon 
substrate. A positive (pointing to silicon substrate) spontaneous polarization 0.26 

C/w? corresponds to a surface potential = 0.962 eV, while P^- ~ —0.26 C/m^ cor¬ 
responds to a surface potential (ps = —0.4346 eV. The depolarization helds through the 
ferroelectric oxide are calculated with the equation (19): 


E, 


r, - / (Q.) 


Qz^z 


(24) 
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FIG. 2. Energy surfaces and their two dimensional projections for MOSFET systems with 
electric oxide BaTiOa of different sizes. Parameters used in the simulation are given in Ref 
(a) dx = 400 nm and dz = 200 nm, polarization in plane favored. Polarizations corresponding to 
the four local minima are marked as Px+, Px-, Pz+ and Pz-] (b) dx = 400 nm and dz = 800 nm, 
polarization out of plane favored. Polarization dynamics is marked with the dashed lines, path 1: 
polarization rotation; path 2: polarization inversion. 


ferro- 

. Q. 


|E.(C,+)l = ;r 

dz 


0.962 + 


0.26A, 


>\Ez{Pz-)\ = ^ 
dz 


0.4346 + 


0.26A, 


(26) 


The depolarization field for positive polarization is larger, and this explains why on the 
energy surface with no gate voltage, a negative polarization is more favorable than a positive 
one. 

Besides, the graphs also demonstrate the known relation that the thicker the ferroelectric 
oxide is in one direction, the more stable the polarization in this direction. As shown in 
equation (19), if the thickness overwhelms the screening length, the potential drop in the 


electrodes can be neglected 


33| . As a result, the electric held through the ferroelectric oxide 
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Charge Density vs. Surface Potential 



FIG. 3. The relationship of surface potential and interface charge density in the doped silicon 
substrate. 


decreases, making the polarization in this direction more favorable. 

These results also illustrate that we can modulate the global minimum by adjusting 
the three-dimensional size of the ferroelectric oxide. An energy surface we are particularly 
interested in possesses the global minimum for P^. When the gate voltage is applied, the local 
minimum corresponding to becomes deeper and polarization rotates to the 2 ; direction. 
After the gate voltage is turned off, the polarization relaxes back along the x direction. 
Meanwhile, the depth of the local minimum for Pz+ is close to that for P^. In such a 
situation, a relatively small applied gate voltage Vg could induce polarization to rotate from 
the X direction to the z direction. The channel current strongly depends on the interface 
charge density, which is approximately equal to the polarization in the z direction. 


.r depends on ^ 

Id, -> Q, ~ P: 


(26) 


Here, we provide guidance about how to select the optimal widths of the ferroelectric 
oxide, in order to make the polarization rotation likely to occur. First, in order to make the 
polarization orient in the x direction without gate voltage, the depolarization field for P^ 
should be smaller than the one that corresponds to P^_, 


Ex 


1 

dx 


2 X 0.26Aa; 


< \Ez{Pz-)\ 


1 


0.4346 + 


0.26A^ 


(27) 
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Typically, the dielectric constant and screening length of the noble metal electrodes are 
£e = 2 and A =0.4 A 261. For these values, we have the criterion 


4 

dx 


< 0.87 


(28) 


In order to have a programmable device, when the applied gate voltage Vg returns 0, the 
polarization should spontaneously return from the 4+ position to the minimum for Px on 
the energy surface. According to the Landau-Khalatnikov equation, the position on the 
energy surface should be a saddle point. 


d^G 


< 0 


Pz=P.+ 


+ ai2P^_f_ + aii2Pz+ + 


2Aj; 

dx^O^e 


< 0 


(29) 


The value of Pz+ increases with thickness in the z direction, since a thinner him means a 
larger depolarization held which suppresses the ferroelectricity. The Pz+ — dz relationship is 
shown FIG. 4. 



FIG. 4. The Pz+ vs. dz plot. Only when the thickness in the z direction dz is below 140 nm, Pz+ 
is smaller than 0.223 C/m^. 


ai + ai2Pf+ + ttii24+ < 0 0 < P^+ < 0.223 C/xa (30) 


Therefore dz < 140 nm is a necessary condition for the polarization rotating back to the x 
direction. 


oil T 


2\x 

dx^O^e 


< 0 ^ 4 > 167 nm 


(31) 
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According to the analysis above, in this study, the BaTiOs dimensions are selected as dx = 
400 nm and dz = 100 nm. 

Hysteresis loops with different values of 70 are calculated and shown in FIG. 5. It demon- 



FIG. 5. Hysteresis loop of the out-of-plane polarization with different values of 70 . 


strates that if 70 is too large, the out-of-plane polarization cannot reduce to 0 and the 
device is not ready for the next program cycle. From our simulation, the threshold 70 for 
out-of-plane polarization returning to 0 completely is around 1.0 x 10^ m/F. 70 is not only 
frequency dependent as shown in equation ( 22 ), but also dependent on the resistance in the 


circuit 34 , l35| , since polarization dynamics is accompanied by screening charge transmission 


26 


36| . Therefore, in order to make 70 in the acceptable range and to have a short switching 


time, the resistance in the circuit should be low. 

To evaluate the performance of MOSFET, The drain-source current Ids and gate voltage 
Vg relationship is calculated based on Pao-Sah double integral. The simulated hysteresis 
loop and Ids — Vg curve for 70 = 1.0 x 10"^ m/F are shown in FIG. 6 (b). 

From the simulation, it can be seen that the on/off ratio of the channel current is large, 
which means that this device is extremely suitable for logic technology. This large on/off 
ratio results from spontaneous polarization rotation, because the spontaneous polarization 
attracts screening charge as free carriers, leading to a large on-current. The segments in 
the Ids — Vg curve circled by dashed lines possess subthreshold swings S lower than 60 
mV/decade. For the segment with Ids and Vg increasing. S' = 53 mV/decade, and the S of 
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FIG. 6 . (a) Hysteresis loop of the in-plane polarization and out-of-plane polarization P^, 

effective polarization dynamic parameter 70 = 1.0 x 10^ m/F; (b) Ids — Vg curve of MOSFET. In 
the circled part, the inverse slope swing is lower than 60 mV/decade. 


the decreasing segment is even lower. This is because, as the gate voltage Vg increases and 
exceeds the threshold voltage, the polarization rotates and boosts free carriers in the silicon 
channel, inducing a steep increase of the channel current. From FIG. 6, we can see that the 
steep change of the channel current is accompanied by polarization reorientation. 

S can be expressed as [17 1 




dV„ 


dV„ 


dps 


d{\og^Q Ids) dpsd{\og^Qlds) 


(32) 


During the polarization reorientation period, the polarization changes suddenly from an in¬ 
plane one corresponding to zero surface potential to a positive out-of-plane polarization, 
which maintains a large surface potential as demonstrated in FIG. 3. The surface potential 
is boosted as 


dps 


< 1 , 


(33) 


causing S to break the 60 mV/decade limit. 

Gompared with polarization inversion, polarization rotation possesses many advantages 
for electronic device applications. First, as shown in FIG. 2, in the polarization rotation 
process^ much lower energy barrier is overcome, leading to a lower polarization rotation 


voltage 


37| |: Secondly, the working state of the MOSFET can be modulated by a unidi¬ 


rectional gate voltage. This simulation and the guidance about designing the ferroelec¬ 
tric oxide size can be also extended to other types of channel, such as quantum well and 


graphene 


38l-l4l|. The only part that needs to be modihed according to the electric 
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properties of new channels is the surface potential-interface charge density relationship 

^s = f (Qz) (34) 

In the simulation, the focus is BaTiOa, but this analytical model can also be applied to 
other ferroelectric oxides, such as PbTiOs and PbZri_ 3 ;Tia, 03 . PbTiOs possesses a larger 


energy barrier in the polarization process compared with BaTiOs 


42l |. Therefore, a larger 


applied gate voltage is needed or we should use PbTiOs with smaller dimensions. Also, 
a single-domain ferroelectric oxide is assumed in this model. However, the effect is not 
limited to a single crystal. When a gate voltage is applied, polarization in the z direction 
increases in different grains and hnally the polarization becomes approximately uniform. 
After the voltage is removed, the polarization relaxes back to the plane. Multiple domains 
may form in each grain, but the polarization distribution in plane has little effect on the 
channel conductance. 


IV. CONCLUSION 

In summary, the polarization distribution in 3D and the electrical properties of the elec¬ 
trodes and the silicon substrate were highlighted in this LGD-theory-based model. Our 
model demonstrated that polarization reorientation can modulate the drain-source current 
effectively. Besides, the choice of electrodes and the dimensions of the ferroelectric oxide are 
key factors in determining the performance of a MOSFET with depolarization helds. With 
proper selection of the thicknesses, held effect transistor with low operating voltage and fast 
switching can be achieved by the polarization reorientation of the ferroelectric oxide. 

ACKNOWLEDGMENTS 

Y. Q. would like to acknowledge the support of the National Science Foundation, under 
grant DMR-1124696. A. M. R. acknowledges the support of the Department of Energy, under 
grant DE-FG02-07ER15920. Both authors thank the National Energy Research Scientihc 
Gomputing Genter for their computational support. 


[1] S. Mathews, R. Ramesh, T. Venkatesan, and J. Benedetto, Science 276, 238 (1997). 


15 




[2] R. C. Naber, C. Tanase, P. W. Blom, G. H. Gelinck, A. W. Marsman, F. J. Touwslager, 
S. Setayesh, and D. M. De Leeuw, Nat. Mater. 4, 243 (2005). 

[3] J. Hoffman, X. Pan, J. W. Reiner, F. J. Walker, J. Han, C. H. Ahn, and T. Ma, Adv. Mater. 
22, 2957 (2010). 

[4] H. Fu and R. E. Cohen, Nature 402, 281 (2000). 

[5] H. Kimura, T. Hanyu, M. Kameyama, Y. Fujimori, T. Nakamura, and H. Takasu, IEEE J. 
Solid-State Circuits 39, 919 (2004). 

[6] L. Liu, V. Narayanan, and S. Datta, Appl. Phys. Lett. 102, 053505 (2013). 

[7] Y. Chen, Y. En, Y. Huang, X. Kong, W. Eang, and X. Zheng, IEEE Electron Device Lett. 
33, 110 (2012). 

[8] H.-P. Chen, V. C. Lee, A. Ohoka, J. Xiang, and Y. Taur, IEEE Trans. Electron Devices 58, 
2401 (2011). 

[9] L. Bellaiche, A. Garcia, and D. Vanderbilt, Phys. Rev. B Rapid Comm. 64, 060103(R) (2001). 

[10] J. Hlinka and P. Marton, Phys. Rev. B. 74, 104104 (2006). 

[11] V. Nagarajan, A. Roytburd, A. Stanishevsky, S. Prasertchoung, T. Zhao, L. Chen, J. Mel- 
ngailis, O. Auciello, and R. Ramesh, Nat. Mat. 2, 43 (2002). 

[12] 1. Batra, P. Wurfel, and B. Silverman, Phys. Rev. B 8, 3257 (1973). 

[13] N. Sai, A. M. Kolpak, and A. M. Rappe, Phys. Rev. B Rapid Comm. 72, 020101(R) (2005). 

[14] W. A. Al-Saidi and A. M. Rappe, Phys. Rev. B 82, 155304 (2010), 

[15] M. A. MendezPolanco, 1. Grinberg, A. M. Kolpak, S. V. Levchenko, C. Pynn, and A. M. 
Rappe, Phys. Rev. B 85, 214107 (2012). 

[16] S. Datta, R. Bijesh, H. Liu, D. Mohata, and V. Narayanan, in Reliability Physics Symposium 
(IRPS), 2013 IEEE International (IEEE, 2013) pp. 6A-3. 

[17] S. Salahuddin and S. Datta, Nano. Lett. 8, 405 (2008). 

[18] C. Krowne, S. Kirchoefer, W. Chang, J. Pond, and L. Alldredge, Nano. Lett. 11, 988 (2011). 

[19] N. A. Pertsev, A. G. Zembilgotov, and A. K. Tagantsev, Phys. Rev. Lett. 80, 1988 (1998). 

[20] V. G. Koukhar, N. A. Pertsev, and R. Waser, Phys. Rev. B 64, 214103 (2001). 

[21] T. Q. Ngo, A. B. Posadas, M. D. McDaniel, C. Hu, J. Bruley, T. Y. Edward, A. A. Demkov, 
and J. G. Ekerdt, Applied Physics Letters 104, 082910 (2014). 

[22] R. R. Mehta, B. D. Silverman, and J. T. Jacobs, J. Appl. Phys. 44, 3379 (1973). 

[23] A. M. Kolpak, N. Sai, and A. M. Rappe, Phys. Rev. B 74, 054112 (2006). 


16 



[24] P. Wurfel and I. P. Batra, Phys. Rev. B 8, 5126 (1973). 

[25] J. Junquera and P. Ghosez, Nature 422, 506 (2003). 

[26] D. J. Kim, J. Y. Jo, Y. S. Kim, Y. J. Chang, J. S. Lee, J.-G. Yoon, T. K. Song, and T. W. 
Noh, Physical review letters 95, 237602 (2005). 

[27] S. M. Sze and K. K. Ng, Physics of semiconductor devices (John Wiley & Sons, 2006). 

[28] Y. Taur, T. H. Ning, et al., Fundamentals of modern VLSI devices, Vol. 2 (Cambridge Uni¬ 
versity Press Cambridge, 1998). 

[29] L. Landau and 1. Khalatnikov, in Dokl. Akad. Nauk SSSR, Vol. 96 (1954) pp. 469-472. 

[30] V. C. Lo, J. Appl. Phys. 94, 3353 (2003). 

[31] W. Zhang and K. Bhattacharya, Acta Mater. 53, 185 (2005). 

[32] H. C. Pao and C.-T. Sah, Solid-State Electronics 9, 927 (1966). 

[33] J. E. Spanier, A. M. Kolpak, J. J. Urban, 1. Grinberg, L. Ouyang, W. S. Yun, A. M. Rappe, 
and H. Park, Nano Lett. 6, 735 (2006). 

[34] S. Sivasubramanian, A. Widom, and Y. Srivastava, Ultrasonics, Perroelectrics and Prequency 
Control, IEEE Transactions on 50, 950 (2003). 

[35] Y. Zhou, H. Chan, C. Lam, and P. G. Shin, J. Appl. Phys. 98, 024111 (2005). 

[36] S. V. Kalinin, C. Y. Johnson, and D. A. Bonnell, J. Appl. Phys. 91, 3816 (2002). 

[37] T. Qi, Y. H. Shin, K. L. Yeh, K. A. Nelson, and A. M. Rappe, Phys. Rev. Lett. 102, 247603 
1 (2009). 

[38] A. Ali, H. Madan, R. Misra, A. Agrawal, P. Schiffer, J. Boos, B. R. Bennett, and S. Datta, 
IEEE Trans. Electron Devices 58, 1397 (2011). 

[39] C. Baeumer, D. Saldana-Greco, J. M. P. Martirez, A. M. Rappe, M. Shim, and L. W. Martin, 
Nat. Comm. 6 (2015). 

[40] X. Hong, A. Posadas, K. Zou, C. H. Ahn, and J. Zhu, Phys. Rev. Lett. 102, 136808 (2009). 

[41] X. Hong, J. Hoffman, A. Posadas, K. Zou, C. H. Ahn, and J. Zhu, Appl. Phys. Lett. 97, 
033114 (2010). 

[42] R. E. Cohen, Perroelectrics 136, 65 (1992). 


17 



